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Abstract 

The Golberger- Treiman discrepancy A^t = 1 — T^^"^ is related to 
. the asymptotic behaviour of the pionic form factor of the nucleon obtained 

00 ! from baryonic QCD sum rules. The result is .015 < A^t < .022. 
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The Goldberger- Treiman relation (GTR) 

rriNQA = IttG^n (1) 



which relates the nucleon mass m^r, the axial- vector coupling constant in decay 
qa, the Trdecay constant /^r and the n—N coupling constant Gttn is one of the most 
^ . remarkable relations of hadronic physics. Explicit chiral symmetry breaking by 

the quark masses leads to small corrections to the GTR, the Goldberger -Treiman 
: discrepancy (GTD) 

^' Aar-l-^ (2) 

which arises from the coupling of the divergence of the axial vector current to 
the = 0~ continuum. The evaluation of Aqt has been addressed recently in 
the framework of baryon chiral perturbation theory f^. On the experimental side 
qa = 1.267 ± .004 and /^r = 92.42 MeV are known to enough precision and most 
of the uncertainty in Aqt results from the uncertainty in G^at. The most recent 
determination of G^at from A^A^, NN and vrA^ data is by the Nijmegen group 

G^N = 13.05 ± .08 which corresponds to Aqt = -014 ± .009 (3) 

Similar results are obtained by the VPI group 0. Larger values are given by 
Bugg and Machleidt ||^ and by Loiseau et al. 0] 

G^N = 13.65 ± .30 which corresponds to A^r = -056 ± .02 (4) 
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The result of theoretical calculations at the loop level do not account even 
for the smaller value given by eq. (|^) in a parameter free way. 

The evaluation of the GTD involves the integral over the imaginary part of 
the form factor Il{q^) which describes the matrix element of the divergence of the 
axial current between two nucleon states 



{Pip') 



d^At \ N{p)) = U{q').U{p'h,U{p) 
p' — p 



(5) 



Access to n(g^) is provided by the study of the three-point function 



d'^xdS exp{-ipx)exp{iqy) (0 | T^^ix)d,A^{y)^^ (0) \ 0) (6) 



where t = p^, d^A'l^ = i{mu + md){u'^^d) expresses the divergence of the axial 



currents in terms of quark fields and 



N 



eijkdi C-f^dj-f^-f^Uk 



(7) 



are the nucleon currents 

The amplitude contains nucleon double and single pole contributions as 
well as a non singular contribution of the continuum 



r{t,q') = {jj) 



X%U{q^)mN 



+ ■ 



(8) 



{t-m%y {t-m%) 

where c is the unknown coefficient of the single pole contribution and Xn repre- 
sents the coupling of the nucleon to its current 

(0 I I P) = XnU^ (9) 

and where we have limited ourselves to the tensor structure for simplicity. 
Any other choice is of course a priori valid provided it leads to stability of the 
calculation as will be shown to be the case here. 

The next step is to evaluate T{t, q^) in QCD. To this end use is made of the 
operator product expansion of the currents entering in eq.(^). The lowest dimen- 
sional operators, which provide the dominant contributions at short distances are 
the unit operator and the operators qq and G^yG'^''^(= GG). As we shall only use 
the coefficient of 75^ in the expansion of the currents entering in the even 
dimensional operators 1 and GG will he multiplied by the small quark mass 
and their contribution will be greatly reduced as compared to the one of the odd 
dimensional operator qq. The contribution of the latter has been evaluated in the 
third of ref. M 



-2mq{qq) 



ln(-t) + lln(-t)-lt/o(t,g2 



+ ■ ■ 



(10) 
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^o(t,?^)= / : dx^^^r -x{l-x)q' ^^^^ 



with 



t — x{l — x)q'^ 

and —2mq{qq) given by the Gell-Mann-Oakes-Renner relation 

-2mg{qq) = f^ml , = .0924GeV (12) 

Expression (|1^) holds for both t and q^ in the deep euclidean region. The 
next step is to extrapolate to the nucleon mass-shell, i.e. to obtain 11'^*-'^ (g^) 
from expressions (|^), ( pUj ) and the analytic properties of r(t,q'^). This extrap- 
olation to the mass shell is done over a large interval of the variables t and q^. 
The method of QCD sum rules provides a tool for such an extrapolation where 
the approximations are well defined and where numerical stability of the result 
provides a useful check of their validity. Note that the smallness of Aqt results 
from the smallness of the quark masses and is not calculated as the difference of 
large numbers. If this were the case the method of QCD sum rules would not be 
accurate enough to be reliable. 

For fixed at a large negative value T{t, q^) has a cut on the positive t-axis 
starting at t = {rriM + m.„Y in addition to the nucleon pole structure exhibited in 
eq. (H). Consider now the Laplace type integral |jTO| 2^ j^dt exp (;^) r(t, q^) in 



the complex t-plane over a closed contour c consisting of a circle of large radius 
and two straight lines above and below the cut which run from threshold to 
i?, is the usual 'Borel mass' parameter. The exponential provides convenient 
damping of the contribution of the integral over the continuum. This contribution 
is of course unknown and provides the main uncertainty in the QCD sum rule 
approach. It could be greatly damped by decreasing the parameter M^, but, as is 
well known this enhances the contribution of the higher order unknown terms in 
the operator product expansion. We hope to obtain an intermediate range of 
for which the contribution of the continuum and that of the higher order terms 
are both negligible. If these approximations are adequate this will show up in 
the stability of the result. We shall find out that this is the case. On the circle 
r is well approximated by r**^'-^°, except possibly for a small region near the real 
axis. 

We obtain then 



27r2A^exp [-^] rriN 



x(l—x) 



(13) 



with 



R 



^^i^)= 1 x'exp{-x)dx (14) 
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n(s = g^) is an analytic function in the complex s-plane except for a simple 
pole at s = and a right hand cut running along the positive real axis from 
s = 9m^ to oo. 

U{s) = -^^-^^f ^ + ■ ■ ■ (15) 
s — ml 

furthermore 

n(s = 0) = 2mNgA (16) 

Consider next the integral J^, ^{s — m''^)Yi{s) where m' is a mass parameter 
and d is a closed contour consisting of a circle of large radius i?'and two straight 
lines above and below the cut which run from threshold to R'^. Cauchy's theorem 
implies 

-m'^ ■ 2mN Qa - 2f^G^N{ml - m'^) 
= i r ^(s-m'2)Imn(s) + — /-(s-m'2)nQC°(s) (17) 

where we have used n(s) = n'^^^(s) on the circle. 

The first term on the r.h.s. of eq. ([T7|) represents an integral over the unknown 
continuum. As m'^ is varied between threshold and R this integral changes sign 
which implies that it vanishes for some value of m'^ which we adopt. Because 
m' is an unknown parameter that we shall vary within reasonable limits it is 
superfluous to include any contribution of the continuum near threshold. The 
GTD then follows from eq.(|l^) 

And when expression ([131) for n**^'-^° is used 



m. 



J dx J dyexp{-y) ( 



4 7o Jo \x{l-x) M2 



+ c"M' 



Xlf is obtained in a similar fashion from a study of the nucleonic two point 
function f d'^x exp{iqx){0 \ T \E'(x)\I'(0) | 0) with the result 

f-ml\ ^ f R\ a^GG, ,,2r. f R\ 32 

(^j = -B. (^j - y(^)^'^^« (iji) + 

(20) 



^R' need of course not be equal to R but they are of the same order and any resonable 
difference between them results only in negligible numerical effects so we take R' = R to 
simplify the notation. 
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The choice of in eq.(]T9|) as well as the consistency of the method is dictated 
by stability considerations. If there are values of small enough to provide 
adequate damping of the continuum and large enough to justify the neglect of 
the contributions of higher order condensates in the operator product expansion 
this should show up in the stability of expression (f^). This means that the first 



term on the r.h.s. of (p!9[) should show a linear behaviour which compensates the 
linear variation of c"M^ in some intermediate range of M ^ (Note that the curve 
need show no horizontal plateau, this happens only if c" = 0). The value of m'^ 
is expected to be close to (albeit smaller because of the weight factor ^) the 
maximum of the 7r'(1.7GeV^) bump. It seems reasonable to vary it in the range 
IGeV^ < m'^ < l.SGeV^. For the gluon condensate we use the standard value 
(^osGG) = .012 GeV^ For {{qqY) the choice (qq)'^ (vacuum saturation hypothesis) 



is usually made but as this seems to be too stringent an assumption |]T0|, we take 
{qq'^) = /5((/g)^. Varying (3 between 1 and 3 has no noticeable effect on the result. 

In the figure the first term on the r. h. s. of eq. ([T9| ) is plotted against for 
m'^ = IGeV^ and P = 1. It clearly exhibits a slow linear variation in the range 
.SGeV^ < M2 < l.SGeV^ which gives 

Act = -022 

varying m'^as discussed above yields finally 

.015 < Act < -022 (21) 

which is consistent with the value given by eq. (^ and clearly favours the smaller 
value of GvrAT- 

It is finally worth investigating the possibility that the value of the quark 
condensate (qq) is much smaller than what results from the GOR relation eq.([l^. 
This is the case for example in "generalized Chiral Perturbation theory" |Tl|. We 
would then have 

2 

Agt^^ or .10<A(,r<.14 (22) 
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